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We analyze the order corrections to the single inclusive jet cross section in lepton-nucleon 
deep inelastic scattering. The full calculation is done analytically, in the small cone approximation, 
. . . obtaining finite NLO partonic level cross sections for these processes. A detailed study of the 

VO ' different underlying partonic reactions is presented focusing in the size of the corrections they get 

' at NLO accuracy, their relative weight, and the residual scale uncertainty they leave in the full 

, cross section depending on the kinematical region explored. The dominant partonic process in very 

forward jet production is found to start at order , being effectively a lowest order estimate, with 
^ ' the consequent large factorization scale uncertainty, and the likelihood of non-negligible corrections 

Q , at the subsequent order in perturbation theory. 
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, Introduction 

oo 

Over the last thirty years, the DGLAP ^ approach to parton dynamics has demonstrated itself as the most adequate 
tool for the description of the energy scale dependence of a variety of lepton-nucleon, and nucleon-nucleon processes 
ly-j , over a wide kinematical range. Surprisingly, not just this approximation, but the lowest order in perturbation theory 
• within this approach (LO) gives fairly accurate estimates for paradigmatic processes such as inclusive deep inelastic 
scattering (DIS), provided an energy or momentum scale of a few GeV characterizes them. The following order (NLO) 
^ I often represents small corrections, required for precise comparisons, but not for the broad picture. 
I " In the few last years high precision DIS experiments, with a wide kinematical coverage, and the ability to measure 
Oh! less inclusive processes, such as those performed by the ZEUS and HI collaborations at HERA, have extended the tests 
on the dynamics of partons to the limits of their kinematical reach, looking for signatures of dynamics complementary 
to that described by the DGLAP approach. Illustrative examples of these tests are the measurements of final state 
hadrons and jets 0,3 produced in DIS processes in the forward region, for which the LO DGLAP description fail 
to reproduce the data by an order of magnitude, and even NLO estimates fall short. 

In a recent analysis Q, it has been shown that the striking failure of the LO description in the case of forward 
Ci I hadrons by no means implies the breakdown of the DGLAP dynamics, but just the inadequacy of the LO picture, 
which simply does not include the dominant contribution to the measured cross section: the process in which an 
initial state gluon is knocked out from the nucleon, and also a gluon fragments into the detected final state hadron. 
Indeed, the NLO approximation, which takes into account these contributions, reproduce nicely the data 0, 0, Q- 

In the case of forward jets, the situation seems to be more compromised, because not only the LO estimates fail, 
but also NLO estimates fall short by a factor of two of the data Q • In reference this feature together with the large 
scale dependence of NLO calculations, has been taken as indicative of the importance of higher order corrections. In 
order to improve the understanding of this situation, in the present paper we compute the order corrections to the 
single inclusive jet cross section in lepton-nucleon deep inelastic scattering. 

We perform the full calculation analytically, factoring out explicitly the remnant initial state singularities into the 
definition of parton densities, and dealing with those of final states in the small cone approximation, which assumes 
that the cone of the jet is narrow This approximation allows us to translate straightforwardly previous results on 
hadroproduction |^ to the case of jets, avoiding a rather cumbersome calculation and delicate numerical treatments 
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for dealing with the coUinear singularities. In Section I wc outline the technical framework required for the small cone 
approximation (SCA) in the case of DIS. 

The small cone technique approximates the full result obtained by NLO Monte Carlo generators, providing fairly 
good estimates for the cross sections with differences typically smaller than the theoretical uncertainties of the NLO 
estimate in the relevant kinematical range. We investigate in section II the range of validity of the small cone 
approximation against the full Monte Carlo NLO parton generators. 

The finite cross sections obtained within the small cone approximation can be implemented in fast and stable 
codes, and simple to use. Exploiting the flexibility of them, in Section III we perform a detailed study of the different 
underlying partonic reactions and their main features: the size of the corrections they get at NLO accuracy, their 
relative weight, and the residual scale uncertainty they leave in the full cross section depending on the kinematical 
region explored. 

Our main conclusion is that, as in the case of hadroproduction, the dominant partonic process in the most forward 
jet region accessed yet is the one with a gluon in the initial state and also a gluon as the main seed of the jet. 
This process starts at order a^, and makes the NLO effectively a lowest order estimate, with the consequent large 
factorization scale uncertainty, and the likelihood of non-negligible corrections at the subsequent order in perturbation. 



I. SINGLE-JET-INCLUSIVE DIS CROSS SECTION IN THE SCA. 



In this section we outline the calculation of the single-jet-inclusive DIS cross section within the small cone approxi- 
mation. This technique has been proposed Q and used in computations of unpolarized single inclusive cross sections 
, and more recently has been extended and validated for polarized proton-proton collisions ^3 • 

The SCA can be thought as an expansion of the partonic cross section in terms of the half aperture S of the cone 
over which the radiation around a given final state parton is integrated. This cone defines the jet at partonic level 
and the integration over it regularizes all final state coUinear singularities. Keeping the most significant terms in 
this expansion, which is given by alog{S) -1-6-1- 0{S^) and neglecting 0{S^) and higher contributions, the expansion 
approximates surprisingly well the full results of the cone algorithm for cone radius up to i? ~ 0.7 in applications for 
proton-proton collisions . 

The main advantage of the approach in the present case is that it is possible to relate analytical results on one- 
particle- inclusive cross sections calculated previously 01, where all coUinear divergences have already been canceled 
or factorized into parton densities and fragmentation functions, to the single-jet-inclusive cross section, for which the 
SCA gives an analytical expression as a function of the cone aperture. As in reference Q, we restrict the discussion 
to the case of non-vanishing transverse momentum. The limit of vanishing transverse momentum requires a more 
involved treatment of coUinear singularities |llL IT^ . Schematically, we start with the finite one-particle-inclusive 
cross section, we undo the factorization of final state coUinear singularities into fragmentation functions, which of 
course is not pertinent in the case of jet cross sections. Then, we add contributions that are not accounted for in the 
one-particle-inclusive fragmentation scheme, i.e. jets formed by two partons. The result is completely finite and for 
phenomenological applications can be convoluted with parton densities just as in the case of the one particle inclusive 
cross section. In practice, the approach is equivalent to having defined effective jet fragmentation functions that 
factorize the final state singularities. 



A. Kinematics 



In order to implement the approach in the case of DIS, we start defining the total cross section for the partonic 
process a + I 1 Jet + I' + X with n partons in the final state: 



dxBdQ^ ~ ^xlSj 



1 



(1) 



where p defines the jet momentum in terms of the momenta of the n partons, ^ is the momentum fraction of the 
parent hadron carried by parton a, Sh = (-P + 0^ the squared energy of the collision in the lepton-proton center 
of mass, and Xb the usual Bjorken variable. The index i = M,L stands for metric and longitudinal contributions 
respectively. Sn is a symmetry factor that account for identical partons in the final state. 

Since the last two factors between square brackets in equation |^ are Lorentz invariant, we can evaluate the latter 
in the hadronic center of mass while the former in partonic center of mass. These frames are defined hy P + q = 
and Pa + Q ^ 0, respectively, with P and pa ~ ^ P the momenta of the initial state proton and parton respectively. In 
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the hadronic center of mass frame, it is convenient to use the transverse energy, E±,and pseudo-rapidity 77, while in 
the partonic frame, we use the same variables as for the one-particle-inclusive case 



y 

Y 



U 



z = 



{Q'^ + s){s + t + u) _ s + t + u 

s{Q^ + s + u) ~ s (1 - y) ' 
(Q2 + S)iS + T + U) _ S + T+U 



S{Q'^ + S + U) 



S{l-Y) 



(2) 
(3) 



with 



S={€I+p^f 



Replacing the jet phase space in eq. ^ we find 



da 



(«) 



S=iq + P)\ 



ttE, 



dxBdQ'^ ^ J dx B dQ'^ dy dz TT ^ {1 — y) 



dE± drj , 



(4) 
(5) 
(6) 



(7) 



where we have introduced 



dxB dQ dy dz 



(8) 



and Cn = I [x\ S\ Sn). We omit from now on the index i, since all the results that follow are valid for both metric 
and longitudinal components. Notice that equation JSJ is valid in d = 4 dimensions. Convolving the partonic cross 
section with appropriate parton densities faiO^ the hadronic cross section can be written as 



1 



1 



1 



2E^dxBdQ^dE^dTj Q^ + Sl-yJo 1 



dz fa{^) do'i'^ij. 



^ dxB dQ"^ dy dz 



At variance with the one-particle-inclusive case, the variable y is completely determined by E±_ and 77, 

y = Y=^e-^ 
V ^ ^1/2 e . 

while z is not fixed. Equation jnj coincides with equation (13) in when Dj{() = S{1 ~ C). 



(9) 



(10) 



B. Jet Contributions 



In addition to the jet analog of the onc-particle-inclusive contribution, where the detected final state jet is originated 
from the fragmentation of just one final state parton, while the other partons play as spectators and are integrated 
over, we have to consider additional contributions in the full jet cross section. The former will be denoted as (Ta^i, 
with initial state parton a and fragmenting parton i. The latter include, on the one hand, contributions accounting 
for the situations in which the jet is formed by two partons which are denoted as (Ta->ij- On the other hand, we have 
to subtract configurations in which the cone that defines the jet contains two of the final state partons, that in the 
one particle inclusive case were nevertheless classified as cfa^i- These contributions will be denoted as as a^^^^^jy 
Discriminated by their partonic content, we have the following ©(a^) contributions: 

• 1^ 9 + 9 + q- 

da^^Sl - 2 do-5_g(g) 

• q ^ q + q' + q': 

+ [da-q^q' - C^0"g^9'(g')] + [daq^q' - C^O-^^q' (g')] + dcFq^qiqi , (12) 



[2daq^g -2daq^gi^g) 



2daq^g(^q)\ 



2 da, 



da, 



g-*aa ' 



(11) 
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[d&g^g - 2 d&q^q^q^] + 2 d&g^qg , (13) 



• .9 ^ .9 + <? + 9: 

[dag^g - dag^g(^g) ^ dOg^gi^^-^ + g ^ g - ^^^^^(g)] + {dOg^g - dOg^gi^g^ + d(7g^gg + dOg^gg . (l4) 

where the &a^i terms are those aheady taken into account in the one-particle- inclusive cross section, with the adequate 
combinatoric prefactors. Notice that cross sections with identical partons in the final state get a symmetry factor, so 
we have to add them to the corresponding matrix elements, before factorization. For instance, all the matrix elements 
involved in the first and third reaction have to be multiplied by 1/2 due to the presence of two gluons or two identical 
quarks in the final state. 

In the small cone approximation it is customary to neglect 0(S^) contributions in the cross section, which means 
we can approximate the matrix elements associated to two partons forming a jet with the corresponding coUinear 
limit. This simplifies drastically the calculation because in this limit the 0(0?^ 1 — > 3 matrix elements factorize into 
0{oLs) 1 — > 2 matrix elements times LO splitting functions. In the list above, we have skipped terms that vanish in 
this approximation. For example, in the second reaction (7 ^ g + + g', the matrix elements vanish in the limit in 
which quark q is coUinear with any of or g'. 

Restricting ourselves to the limit when partons j and k become coUinear, for both daa_tj(k) ^^nd d&a^jk, the 
corresponding matrix elements are given by ^3 (d = 4 + e): 

\Ma^ijk\^ ^ —4nfi-'a, \Ma^^J\'' Pj^Ma,^). (15) 

Sjk 

where the momenta of the parent parton, J, is defined in terms of those of partons j and k: 

pi^^ap^j + k^^-^-^, fcM_ JL^^^, (16) 

a 2pj ■ n 1 — a 2pj ■ n 

with pj ■ k± — n ■ k± — and p^ — — 0. Consequently 

PJ = ^-^^^^^^ - k.. (17) 
1 — 2 a 

Notice that due to the singularities (poles in e) arising from the sj^ factor in the denominator when integrating over 
the phase space, we need the splitting functions in l|15(l up to order e. The poles in a = 1 in the diagonal kernels 
(J = j) will give rise poles in z = 0. 



C. daa^j(k) contributions 



As we mentioned, these contributions correspond to configurations where the jet is formed solely by parton j, but 
parton k lies within the cone. They can be computed along the lines of what has been done in reference [lol | but now 
for DIS kinematics. In the partonic center of mass frame wc have 



Sjk = 2pj -pk = 2Ej Ek (1 - coaO-jk) , 



,1/2 



(i-^(i-y)), 



EkC^—z{l-y), 



while the coUinear limit implies a = \ — z{l — y). Then, we can write the matrix element for daa^j^k) 



ijk\ 



\\kl-zil-y) 



1 



zil-y) E^{l-cosejk) 
For the phase space we have. 



P.j^jk{a,e) An "as 



Ma~*,j {s,y' = 



1- zil-y) 



dPS 



(3) 



(4 7r)4 ^^^) 



(477) 



d-4 



9,k dcos 9,k{l -yfz{{l- yf y z^ (1 - z)) 



(d-4)/2 



(18) 



(19) 



dydz, (20) 
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and combining with the matrix elements in H19() . we have 



da, 



a~*j{k) 



1 



dx B dQ^ dy dz Sn \ 2 



(21) 



where 



a— — 



x|5|(4^)2 e r2(l + |) l-z{l-y) 



X ^ 



Ma-.ij [s,y' = 



l-z(l-y) 



(22) 



Notice that with this notation, the variables y and z in !Fa^ij{s^ y, z) are defined in terms of the second parton in the 
final state J. As mentioned above, the poles in a = 1 in the diagonal kernels, appear as poles in z = (notice that 
?/ = 1 is protected by the l — y factor coming from the phase space) . This poles have to be prescribed as usual, leading 
to double poles in e, together with the appearance of 6{z) and 'plus' distributions in z = 0. The double poles cancel 
with terms coming from daa—>jk, whereas the remaining 1/e poles cancel terms generated when the fragmentation 
function factorization of the one-particle-inclusive cross section is undone. 



D. daa~,jk contributions 

These terms account for jets formed by two partons, j and k. In the coUinear limit (|15|l we have: 



a L — a 



(23) 



where a is again the momentum fraction in the splitting function. The phase space in this case is 



dPS^^^ 



(4 7r)3 T^^) 



W '5(z)(y(l-y))('^-4)/2 (^^y ' ^a''-^ da su,'^-' 9, de.dydz, 



(24) 



where we have chosen the frame in such a way that the jet momentum is oriented in the z axis. This result exhibits 
some differences with the previous case. Due to the presence of the S{z) function, the daa~*jk pieces only contribute 
at z = 0. On the other hand, now we have to perform explicitly the integration over the argument of the splitting 
functions, a. In addition, we see that the only dependence on the angle in the matrix element comes, again, form the 
Sjk denominator. However, now we are integrating over 9j and not over 9jk. As in [lol |. we find 



a -a) 



(25) 



and the integration limit over 9j is given by 6 if Ej < Ek and by (l — a)/aS if Ej > Ek- Performing the angular 
integrals in the partonic center of mass frame, 



da, 



a^jk 



1 



dx B dQ^ dy dz Sn 



Ta^tjS{z) 



da 



e 



a I a""' + Q { a 



(l-aY 



Pj^,fc(a,e), (26) 



E. Cancellation of singularities. 

The last step consists in reverting the factorization of final state collinear singularities already done in the finite 
one-particle-inclusive cross section. The simplest way to implement this step is simply to add to the already finite 
one-particle-inclusive cross section the terms we customary factorize into fragmentation functions, with the opposite 
sign. This provides a check of the whole calculation, as the simple poles still appearing in the correction terms have 
to cancel when we add these contributions. 

The terms that have to be added can be read from the following factorization prescription formula for the partonic 
cross sections at first order, when rcnormalizing fragmentation functions: 



dxB dQ^ dy 



(27) 
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The splitting kernels in H27|l are the regularized ones in 4 dimensions. The Born level cross sections are, in turn, given 
by: 



dxBdQ^dy x% Sj^ 8tt r{l + e/2) \4tt^ 



=/2 



(28) 



The additional factor of 1/^ in the hadronic cross section is simply re-expressed in terms of the new variables and 
factors out everywhere. 

Taking into account all the first order cross sections, we have the following contributions coming from reverting 
final state factorization 



dact '^'d&q^^qy) + ddq^(gq) + d&q^(gg) + ^ (d(T,_ + d& q_,^^g, g,)) 

q' 

+ d(Tg^[qg) + d&g^^gg) + dCTg^^^gg) + dCTg^^gg) , 



where we have introduced the notation da, 



da, 



dx B dQ'^ dy dz 



r^(i + £/2) fM'l) 
r(i + e) 



a~>(jk) ■ 
2 \ 



J^a^U X (1 - y) (1 - z (1 - y))-' P.J-,A^ ~z{l-y)). 



(29) 



(30) 



Notice that we have already added a minus sign to these contributions, as they must be subtracted to recover the 
original, un-factorized, one-particlc-inclusive cross section. We also took into account the fact that at first order, the 
sum in eq. H28|l contains only one term. 

As mentioned, doubles poles must cancel explicitly between the aa^j(k) and aa^jk contributions. The remaining 
single poles must cancel with daa^ijk). To make this cancellations more transparent, it is convenient to group the 
different contributions in eqs. 1)1111 - 1)14^ and the ones in H29II . Starting with the quark initiated reactions and omitting 
daa^j contributions (which are already finite) we have 



2dag^g^g) ~ d^g^^gg) 

2daq^q(g) — 2d(7g-,qg 



-Uf 

— da. 



do'q-,q'{g') + daq_,q,^q,) — dcFq^^qiq,-^ — d(7g_^(g/g') 



2daq^g(^g) -daq^qq -daq^(qq) — U f d&q^ 



q^(.gg) 



For the gluon initiated reactions, we have: 



'^^g^giq) - '^^g^(gq) 



da. 



g~*g{q) 



da. 



g-*{gq) 



da. 



g^q{g) 



da. 



g~*qg 



da. 



'q^q'q' 



(31) 



g^iqg) 



d^g^q(g) ^ d&g^gg - ddg^^gg) 



(32) 



In Appendix I, we list the results corresponding to the cancellations for each of the terms in square brackets, which 
completes the calculation. 



II. COMPARISON WITH MONTE CARLO RESULTS. 



Having obtained the finite expressions for the NLO corrections in the SCA, in this section we investigate the 
accuracy of the approximation. In Figure^we compare the outcome of the SCA for the single jet inclusive DIS cross 
section with a full Monte Carlo NLO calculation of [3. The jets are reconstructed in the Breit frame and the rates 
between both results are computed in the typical kinematic range of forward jet DIS experiments at HERA. Jets are 
defined using the inclusive fc^ cluster algorithm in the Monte Carlo, and a cone radius of i? = 0.7 for the SCA 
value for which the agreement between both jet definitions is maximized. 

The rates are presented as a function of the transverse momentum Et and rapidity -q of the jet, both measured 
in the laboratory frame, and Bjorken momentum fraction xb In both cases we use the MRST02 NLO set of parton 
densities j3 ^-nd we compute at NLO fixing Aqcd as in the MRST analysis so as{Mz) = 0.1197. The rapidity 
variable varies between 2.0 and 3.5, the Bjorken variable spans the interval between 0.0005 and 0.01, the transverse 
momentum of the jet starts at 5GeV, and r ang es from 20 to 100 GeV^. 

As it has been observed in references and |lOl |. for hadronic collisions, even for rather large cone radius the SCA 
gives acceptable approximations within less than a ten percent of the full Monte Carlo result. This is also the case 
for DIS and the rates show a very mild dependence in the kinematical variables. 



7 




FIG. 1: Ratio of NLO SCA estimate and the full Monte Carlo prediction with kr jet reconstruction. 



Certainly, the accuracy of SCA is the better for smaller cone radius, however the error introduced by the approx- 
imation with R ~ 0.7, which is of the order of a 10 %, always underestimating the cross section and with a very 
mild dependence on the relevant variables, is comparable or smaller than the theoretical uncertainty coming from 
the particular choice of the factorization and renormalization scales, characteristic of the NLO corrections in this 
kinematical region. As we show in the following section, this means that we can safely use the SCA results as an 
estimate of the size and behavior of the NLO corrections. For comparison, in the Figure we also plot as a band 
the uncertainty resulting from varying the factorization and renormalization scale by a factor of two in the full Monte 
Carlo result. 



III. PHENOMENOLOGICAL CONSEQUENCES 



Having established our level of confidence in the SCA results, we proceed analyzing the distinctive features of 
the NLO corrections in the forward region. We do the analysis in the typical kinematic region tested by DESY 
experiments, where large higher order effects have been observed. 

The most striking of these features is the size of the NLO corrections as the rapidity of the jets increases. In Figure 
12 we show both LO an NLO partonic level expectations coming from the SCA approach in three different regions 
of rapidity for 5 < < lOOGeV^ and 0.2 < yei < 0.6 as a function of transverse jet momentum. Clearly, NLO 
corrections, which are moderate for central rapidities, become significantly large in the forward region. This feature is 
due to a suppression of the LO contributions rather than to an increase in rapidity of the cross sections. K-factors can 
exceed an order of magnitude there, invalidatin g th e lowest order approximation. We have included for reference the 
data obtained by HI in that kinematical range |l7l | , although one should keep in mind that for a precise comparison 
hadronization effects should also be taken into account. 

Another interesting feature of NLO corrections to be taken into account is the rather large uncertainty these 
corrections show associated with the choice for the factorization and renormalization scales. In Figure [21 we adopted 
^2 — q2 ^j^g factorization and renormalization scale. The choice for these scales is in principle arbitrary; the 
differences found in any perturbative estimate coming from some particular choice for the scale or other, become 
smaller as more terms in the perturbation series are included. In inclusive DIS is the typical choice, while is 
the one favored in jet physics. In Figurc|2|wc plot the uncertainty bands corresponding to vary the scale from to 
E^, the two main scales of the process under consideration. 

With such a large scale uncertainty, any particular choice will probably lead to miss the data at some point. One 
possible choice in these cases is to take the average between them, which leads to an intermediate estimate within the 
band. Again we can see that in the most forward bin, the uncertainty associated to the choice of the scale becomes 
more prominent and may be as large as a factor of two. 

In order to understand the correlation between higher order corrections and rapidity, it is useful to discriminate 
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FIG. 2: Size of NLO corrections as a function of Et for different rapidity regions. 



Ji^ 10- 

Oh 10^ 



-l<x\<0.5 



0.5<r[<1.5 



1.5<x\<2.8 



I 



10 



10 



10 



HI 

NLO elO' , E^] 



10 



10 



10 



Ej. [GeV] 



FIG. 3: Scale uncertainty in NLO corrections as a function of Et for different rapidity regions. 



the different partonic contributions, classifying them as in the one particle inclusive case in terms of the initial state 
parton, i, and the one taken as the seed for the jet, j, as in cTi-^j. The additional contributions coming from the first 
two terms in Eq.(32) were added to the (Tg_yg contribution, the remaining terms in that equation were associated to 
Cg-tq. Contributions in Eq. (31) were taken together with Oq^g. 

In Figure 0] we show the different partonic contributions in the three rapidity regions. While in the central region 
the LO contribution is very close to the full NLO estimate, in the forward region it is significantly smaller. In the 
former region, the cross section is dominated by the Cg^q contributions (an initial state gluon with a quark originating 
the jet), which are already present at LO, while in the latter the dominants are Cg^g, which are pure NLO, shown in 
Figure These contributions start at order aj. so the NLO result is its lowest order estimate. The reason for their 
dominance over the LO is just the kincmatical region chosen which suppress LO configurations. The dominance of 
the (7g^g over cfg^q NLO contributions can be traced back to the negative 'plus' contributions which are proportional 
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FIG. 4: Partonic contributions as a function of Et for different rapidity regions. 



C|n in the case of the former and to CfCa for the latter. 









FIG. 5: Gluon initiated contributions at 0{o?a) 



Since the dominant partonic process in the forward region is accounted, at order a^, only by only its lowest order 
contribution, it is effectively a LO estimate and most probably receives significant higher order corrections. The first 
order corrections for other partonic processes in this kinematic region rise typically to 50% effects, so it would not be 
surprising that the NLO estimate falls short of the data, specially if a more stringent kinematic range is explored. 
This is precisely what ZEUS has reported in preliminary analyses of measurements in the very forward region |l8l | . 

In Figure ini we plot the NLO estimates for the cross section as distributions in different variables together with 
ZEUS preliminary data . The estimate correspond to rapidities between 2.0 and 3.5, the Bjorken variable in the 
interval 0.0004 and 0.005, the transverse momentum of the jet starting at 5 GeV, and the virtuality of the photon 
range from 20 to 100 GeV^. The NLO estimate falls short of the preliminary data, and only allowing a rather large 
scale uncertainty it may be considered consistent with the measurements, specially at small xb- 

Further insight is obtained analyzing the different partonic contributions as a function of rj and xb- In Figure [7| it 
can be noticed that the cTg^g contributions dominate the cross section, specially at low xb where the gluon parton 
density grows dramatically and in the middle of the rapidity range. In these two regions one can expect the first order 
corrections to these processes, starting at NNLO, to be significant. In fact, it is there where the NLO estimate can be 
more distant to the data with a particular choice for the scale, as can be seen when comparing with the preliminary 
data in Ref. |l8l |. At larger xb and rj, (Jg^g contributions decrease, even below the LO contribution, but the other 

(ii) 

NLO contributions keep the total NLO effect very large. Of these contributions, the most prominent is Ug-Zg, that 
corresponds to diagrams, like the ones in Figure |H1 where the two quarks lines have different flavours, but where the 
quark that initiates the jet has the same flavour as the one in the initial state. 
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FIG. 6: NLO estimates against ZEUS very forward preliminary data |I| 




FIG. 7: Partonic contributions in Zeus very forward measurement [l3 | 



Going to even higher rapidities, these last contributions eventually dominate the cross section, with the LO estimate 
being completely suppressed as shown in FigurcEl In this region, with the two dominant contributions being computed 
at the lowest order, the scale uncertainty of the NLO estimate is twice as large as that found in the rapidity region 
of Figure 13 
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FIG. 8: Typical quark initiated contributions at 0{a^) 




3.5 4 4.5 5 0.0004 0.0027 0.005 



FIG. 9: Partonic contributions in an extremely forward region 
IV. CONCLUSIONS 

We have computed the single jet inclusive deep inelastic scattering cross section at C(a^) in the small cone approx- 
imation. We found that this approach approximates the full NLO Monte Carlo results within a 10% accuracy, error 
which is fairly moderate compared to the main source of theoretical uncertainty, the scale dependence. 

As in the case of hadroproduction in deep inelastic scattering, we have found that the dominant partonic processes 
in very forward jet production start at order a^, being effectively a lowest order estimate. As in any lowest order 
calculation, there is a large factorization scale uncertainty which can not be neglected, and it is likely that there will 
be large corrections at the subsequent order in perturbation. Although taking into account this large dependence on 
the choice for the scale, one can bring agreement between data and NLO estimates, the difference between them for 
a particular choice is maximal precisely where the partonic contributions computed for the first time are dominant. 
This feature is expected to be even more apparent at higher rapidities, and the corresponding measurements will 
constitute an obligatory benchmark for the study of QCD at NNLO. 
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VI. APPENDIX 



Here we list the finite results obtained for each of the terms in square brackets in eqs. (31) and (32). 



r(l + e/2) fMl^'^^ 



r(l + e) V s 



T,,2 = -2n/.F,_ 



(1 - y) (1 - z (1 - y))-^ P,^,,(l - z (1 - y)) 
(l~y)i+'z^P,^,,(l-z(l-2/)) 



Tq^3 ~ ^q^gq 



(l-2/)i+^z^P,_,g(l-z(l~y)) 



5{z) I daG{a,e) Pq^qg{a,e) 



r(l + e/2) (Ml^''^ 



r(l + e) V s 



(1 - 2/) (1 - z (1 - y))-^ P,_,,(l - z (1 - y)) 



Tq^4 — —J^c 



q^qg i \ 2 
5(z) 



(l-y)i+^z^P,^,,(l-z(l-y)) 



daG{a,e) [Pg^gg{a,e) + 2nf Pg^qq{a,e) 



r(l + 6/2) fAP^^'/' 
Til + e) 



(1 - y) (1 - z (1 - y))-^ P„(l - z (1 - y)) 



TgS - 



il-yy+^z^Pq^gqil~z{l-y)) 



r(l + e/2) fM^ 



r(l + e) V s 



2 \ «/2 



(1 - y) (1 - 2 (1 - y))- P,^,,(l - 2 (1 - y)) 



^ff,2 — ^g^qq 



{l-yY+^z^Pq-,gq{l-z{l-y)) 



Tirrr (^)'^' (1 - y) (1 - ^ (1 - 2/))"^ ^^,^.,(1 y)) 



— -Fg^qq 



{\-yf^^z^Pq^qg(\-z(\-y)) 



b{z) \ daG{a,e) Pq^qg{a,e) 



r(l + e/2) (MlX^^ 



r(i + e) V s 



(1 - y) (1 - z (1 - y))-^ P,_,g(l - z (1 - y)) 



(33) 



(34) 



(35) 



(36) 



(37) 



(38) 



(39) 



13 



TgA — I 



Jo 
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